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ABSTRACT 


The  active  control  of  acoustic  pressure  in  a, '2^  cavity  with  a  flexible  boundary  (a  beam) 
is  considered.  Specifically,  this  control  is  implemented  via  piezoceramic  patches  on  the  beam 
which  produce  pure  bending  moments.  The  incorporation  of  the  feedback  control  in  this 
manner  leads  to  a  system  with  an  unbounded  input  term.  Approximation  methods  in  the 
context  o||4rp^QIi[)  state  space  formulation  are  discussed  and  numerical  results  demonstrating 
the  effectiveness  of  this  approach  in  computing  feedback  controls  for  noise  reduction  are 
presented.  _ 
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In  recent  years,  the  development  of  new  fuel  efficient  turboprop  engines  has  motivated  the 
development  of  a  comprehensive  active  control  methodology  for  interior  pressure  field  cham¬ 
bers.  The  active  control  of  noise  in  this  setting  has  been  studied  both  in  a  frequency  domain 
setting  [14,  19]  and  from  an  infinite  dimensional  state  space  time  domain  approach  (PDE 
approach)  [2,  6,  7, 12]  with  techniques  often  centering  around  the  generation  of  an  appropri¬ 
ate  secondary  pressure  wave  which  optimally  interferes  with  the  offending  primary  pressure 
wave.  Here  however,  we  consider  a  time  domain  state  space  formulation  in  which  the  active 
control  is  implemented  via  piezoceramic  patches  which  are  imbedded  in  the  boundary  of  the 
acoustic  cavity. 

The  example  we  consider  consists  of  an  exterior  noise  source  which  is  separated  from 
an  interior  chcimber  by  an  active  wall  or  plate.  This  plate  transmits  noise  or  vibrations 
from  the  exterior  field  to  the  interior  cavity  via  fluid/structure  interactions  thus  leading  to 
the  formulation  of  a  system  of  partial  differential  equations  consisting  of  an  acoustic  wave 
equation  coupled  with  elasticity  equations  for  the  plate.  The  control  is  implemented  in  the 
example  via  piezoceramic  patches  on  the  plate  which  are  excited  in  a  manner  so  as  to  produce 
pure  bending  moments.  It  should  be  noted  that  the  incorporation  of  the  feedback  control 
in  this  manner  leads  to  a  system  with  an  unbounded  input  term.  Experiments  are  being 
designed  and  carried  out  at  NASA  Langley  Research  Center  in  which  the  interior  cavity  is 
taken  to  be  cylindrical  with  a  circular  active  plate  and  sectorial  patches. 

As  a  first  step  toward  developing  an  effective  linear  quadratic  regulator  (LQR)  state  space 
control  methodology  for  near  field  acoustic  problems  of  this  type,  it  is  useful  to  consider  a 
simplified  but  typical  model  consisting  of  a  2-D  interior  cavity  with  an  active  beam  at  one 
end  (see  Figure  1).  Here  T  represents  a  perturbing  force  on  the  beam  due  to  an  exterior 
noise  source.  This  in  turn  causes  fluctuations  in  the  interior  acoustic  pressure  field  and  hence 
unwanted  noise.  The  goal  in  the  control  problem  is  to  optimally  reduce  the  interior  pressure 
deviations  by  effecting  a  force  distribution  on  the  beam  that  decouples  the  cavity  acoustic 
response. 

In  Section  2,  a  model  set  of  differential  equations  for  the  problem  is  given  and  the  math¬ 
ematical  framework  needed  to  pose  the  control  system  in  an  abstract  Cauchy  formulation  is 
presented.  Section  3  contains  a  brief  discussion  of  the  theory  of  finite  and  infinite  dimensional 
periodic  optimal  control  problems  while  Section  4  is  devoted  to  the  general  finite  dimensional 
approximation  of  the  control  problem.  Specific  approximation  schemes  are  discussed  in  the 
fifth  section  and  examples  demonstrating  the  viability  of  the  method  are  presented. 
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Figure  1.  Acoustic  chamber  with  piezoceramic  patches. 


2  Mathematical  Model 

When  describing  acoustic  wave  motion  in  a  fluid,  it  is  useful  to  introduce  a  velocity  potential 
<f>  which  is  a  complex- valued  function  satisfying  v{t,  x,  y)  =  —  x,  y)  where  v  denotes  the 
fluid’s  velocity  [15,  16].  If  the  equilibrium  density  of  the  fluid  is  given  by  p/,  the  acoustic 
pressure  p  (the  deviation  from  the  mean  pressure  at  equilibrium)  is  related  to  this  velocity 
potential  by  p(t,x,y)  =  For  acoustic  waves  with  small  amplitude,  both  the 

potential  and  the  pressure  satisfy  the  undamped  first  order  wave  equation  with  uniform 
speed  of  sound  c  in  the  fluid;  hence 

<f>tt  =  c^^(f>  (x,y)  6  n(t)  ,<>  0  . 

The  boundaries  on  three  sides  of  the  variable  cavity  are  taken  to  be  “hard”  walls  thus 
leading  to  the  zero  normal  velocity  boundary  conditions 


V^*n  =  0  (x,i/)€r,t>0 


where  h  is  the  outer  normal.  It  is  assumed  that  the  perturbable  boundary  consists  of  an 
impenetrable  fixed-end  Euler-Bernoulli  beam  with  Kelvin- Voigt  damping.  If  w{i,x)  is  used 
to  denote  the  transverse  displacement  of  the  beam  with  linear  mass  density  pt,,  the  equations 
of  motion  are 


Pb'Wtt  +  =  -pf(f>i(t,x,w{t,x))  -f  f{t,x) 

w{t,0)  = -^{t,0)  =  w{t,a)  =  ■^{t,a)  =  0  t>0  , 


(2.1) 


where  M(<,  i)  is  the  internal  moment  and  /  is  the  external  applied  force  due  to  pressure  from 
the  exterior  noise  field.  For  an  uncontrolled  beam  with  Kelvin- Voigt  damping,  the  moment 
contains  both  strain  and  strain  rate  components  and  is  given  by 


M{t,x)  El  ■ 

The  final  coupling  equation  is  the  continuity  .of  velocity  condition 


=  V(^(t,a:,to(t,a:))  •  n  ,  0<x<a,t>0  (2.2) 

which  results  from  the  assumption  that  the  beam  is  impenetrable  to  fluid.  Under  an  eis- 
sumption  of  small  displacements  (t«(<,  x)  =  w{t, «)  -f-  5  where  to  =  0)  which  is  inherent  in 
the  Euler-Bernoulli  formulation,  the  beam  equation  in  (2.1)  can  be  approximated  by 

Pb'^tt  +  (^>  ^)  =  -‘PAM^^  2:,  0)  +  (j>ty{t,  X,  0)u;]  +  f{t,  x) 

while  (2.2)  can  be  approximated  by 


Wt{i,x)  =  V^6(<,i,0)  •  n  +  {V<l>y{t,x,0)w)  •  n  . 


To  first  order,  these  last  two  equations  can  be  approximated  by  dropping  the  higher  order 
terms  —pf(i>ty{t,x,0)w  and  (V^y(t,x,0)w)  -n.  Then  upon  approximating  the  domain  Q(t) 
by  the  fixed  domain  fl  =  [0,  a]  x  [0,£],  we  obtain  the  approximate  uncontrolled  model 


(x,^)  €  Q  ,t  >  0  , 

V^-n  =  0  (®,j/)€r,<>0, 

~{t,x,0)  = —wt{i,x)  0<x<a,t>0, 
ay 

w(<,0)  = -^(^,0)  =  to(<,a)  = -^(^a)  =  0  t>0  , 


0  <  a:  <  a  , 
t  >0  , 


(2.3) 


^(0,  X,  y)  =  (i>o{x,  y)  ,  w{0,  x)  =  7i;o(x) 

<^t(0,X,J/)  =  (^l(x,t,')  ,  tOt(0,x)  =  tOi(.T)  . 
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For  control  of  structural  vibrations  and  the  acoustic  pressure  field  in  this  model,  s  piezo¬ 
ceramic  patches  are  attached  to  the  beam  as  shown  in  Figure  1.  These  patches  are  excited 
in  a  manner  so  os  to  produce  pure  bending  moments  ([8,  9,  11])  (see  Figure  2).  If  H  is  used 
to  denote  the  Heaviside  function,  the  model  for  the  controlled  beam  can  be  written  as 


5^  (  d'^w  d^w  ^ 
^  dj!p  ) 


+  pj<pt{t,x,0) 


dx^ 


^  Ui{t)  [H{x  -  an)  -  H{x  -  a, •2)])  -f  /(t,  x)  . 


(2.4) 


Here  «,(t)  is  the  voltage  applied  to  the  patch,  is  a  parameter  which  depends  on  the 
geometry  and  piezoceramic  material  properties,  T  is  the  patch  thickness  and  k  is  a.  material 
constant  (see  [8,  9]).  It  should  be  noted  that  the  incorporation  of  (2.4)  into  (2.3)  leads  to  a 
system  with  an  unbounded  input  term  since  it  involves  the  second  derivative  of  the  Heaviside 
function. 


-I- 


Figure  2.  Piezoceramic  patch  excitation. 

To  formulate  this  problem  in  the  context  of  existing  infinite  dimensional  control  theoretic 
results,  it  is  advantageous  to  pose  the  control  system  in  an  abstract  Cauchy  formulation.  To 
accomplish  this,  the  state  is  taken  to  be  ^  =  {<f>,w)  in  the  Hilbert  space  H  L^{n)  x  L'^iFo) 
with  the  energy  inner  product 

((“■)’(  0)„  "  +  ■ 

Here  L^{Sl)  is  the  quotient  space  of  -ver  the  constant  functions.  We  also  define  the 
Hilbert  space  V  =  fP{Q)  x  //^(Fo)  where  /7*(fi)  is  the  quotient  space  of  IP  over  the 
constant  functions  and  ’’^(Po)  =  {^  E  lP{To)  •  ip{x)  =  ip'{x)  =  0  at  x  =  0,  a}.  The  V  inne^ 
product  is  taken  as  (heio  and  below  we  u.se  the  notation  T)  =  ^) 
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Following  the  ideas  used  in  the  theoretical  results  in  [3,  4],  we  consider  the  Gelfand  triple 
V  ^  H  ^V*  with  pivot  space  H  and  define  sesquilinear  forms  cr,:VxV'— >€,  i  =  l,2 

by 

<ri($,  ){')=  f  pfV^-V^du)  +  f  ETD^wD^'qd'y  , 

Jn  Jvo 

(72(^5  ^)  =  /  {cdID^wD^t}  +  pf{<j>r]  —  w()}d'y 
JVo 

where  $  =  {<f),w)  and  ^  =  (4,j/)  are  in  V.  It  can  be  easily  argued  that  the  sesquilinear 
forms  satisfy  the  continuity  and  coercivity  conditions 

lcri($,^)|  <  cal^lvl^lv  , 

Re<r2($,$)  >  =  caliaj^Kp^)  , 

k2($,^)l<C4l$lv|^k 

(for  detailed  arguments  in  a  similar  setting,  see  [1]).  The  control  operator  B  G  C{U^V*)  is 
defined  by 

r  K^k  * 

•'I'o  ■'  t=l 

for  '3?  €  V,  where  Hij{x)  =  H{x  —  Qij),i  =  l,2,"',s,  j  =  1,2  and  {-,-)v,v  ^be  usual 
duality  pairing. 

Finally,  for  F  =  (0,  f/pb)  we  can  write  the  control  system  in  weak  or  variational  form 

Mt),  'I')v.v  +  +  ari{z{t),  '9)  =  {Bu{t)  +  F,  '3?)v.,v  (2-5) 

for  '5  in  V.  The  state  is  given  by  z{t)  =  (^(t,  inV^H.  Since  <7i  and  cr2  are 

bounded,  we  can  define  operators  A2  €  £(V,  V"*)  by 

for  i  1,2.  This  then  yields  the  system 

+  A2Zt(i)  +  Aiz{t)  =  Bu{t)  +  F 

in  V\ 

Continuing  with  our  abstract  formulation,  we  next  write  the  system  in  first  order  form. 
To  accomplish  this,  define  the  product  spaces  V  =  V'  x  and  H  =  V  x  II  with  the  norms 

and 

For  X  =  (^> '!')  ond  0  =  (T,  A),  the  sesquilinear  form  <7  :  V  x  V  — »  €  is  then  defined  by 

a((T,  A),  ($,  'P))  =  -  (A,  ❖}v  +  <T, (T,  +  (TjC A,  .  (2.6) 
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Since  the  duality  product  (•,  •)y,  y  is  the  unique  extension  by  continuity  of  the  scalar  product 
H  xV  to  V*  xV,  it  follows  that  for  appropriate  restrictions  on  0  we  can  write 

=  A),  ($,$))  =  ~(A,$)v.  +  (AiT,$)^.,^  +  (A2A,^)v...v 

=  -(A,$)^  +  (^T  +  .42A,^);, 

=  ((-A,AxT  +  /l2A),(<&,^))« 

=  (->!€),  X)w  - 

The  operator  A:'H—*'His  given  by 


where  dom  =  {0  =  (T,  A)  eH  :  A  €  V,AiT  +  >l2A  G  //},  Ai  and  A2  are  the  operators 
defined  by  <ti  and  <72,  respectively,  and  the  above  calculations  hold  for  0  G  dom  A  (see  [1] 
for  further  examples  concerning  the  definitions  of  operators  and  domains  in  this  manner). 

To  write  the  first  order  system  in  weak  or  variational  form,  let  Z{t)  =  {z{t),zt{t)), 
^(t)  =  (0,F(t)),  and  Su(t)  =  (0,Bu(t)).  The  weak  form  of  the  system  is  then 

+  =  +  (2-8) 
for  X  €  V.  Formally,  this  is  equivalent  to  the  system 

Zt(t)  =  AZ(t)  +  Bu(t)  +  F(t)  (2.9) 

in  W  where  A  is  given  in  (2.7). 


3  Periodic  Control  Problems 

As  noted  in  the  introduction,  ot  *  '•ontrol  problem  is  motivated  by  the  desire  to  reduce 
cavity  pressure  fluctuations  result:,  v  from  the  perturbing  noise  T.  In  many  applications, 
it  is  reasonable  to  assume  that  F  periodic  with  period  r;  hence  an  important  problem 
of  interest  (e.g.,  see  [6])  for  the  system  (2.9)  is  an  LQR  problem  for  f  periodic  disturbing 
force  y.  This  c.'.n  be  formulated  as  the  problem  of  finding  u  G  >i/‘*(0,  t;  f/)  which  minimizes 
a  quadratic  cost  function?’  of  the  form 

m  =  5  l{aZ(t),Z{t))n  +  {Ru{t),u(t))u)  dt 

subject  to  (2.9)  with  Z(0)  =  Z{t).  Since  Z  =  (^,t«,  At,  wt)^>  l-hc  operator  Q  can  be  chosen 
so  as  to  emphasize  the  minimization  of  particular  state  variables  as  well  as  to  create  windows 
that  can  be  used  to  decrea.se  state  variations  of  certai;'  rrequencies.  The  control  space  V  is 
taken  to  be  IR*  if  $  patches  arc  used  in  the  model,  and  is  assumed  that  the  operator  R  is 
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an  s  X  5  diagonal  matrix  where  r,-,-  >  0,  i  =  1,  •  •  • ,  5  is  the  weight  on  the  controlling  voltage 
into  the  patch.  In  the  Ccise  that  B  is  bounded  on  H,  a  complete  feedback  theory  for 
this  problem  can  be  given  as  discussed  in  [10].  Under  usual  stabilizability  and  detectability 
assumptions  on  the  system  as  well  eis  standard  assumptions  on  Q,  the  optimal  control  is 
given  by 

u{t)  =  -R-^B*[YiZ{i)  -  r(t)] 

where  11  is  the  unique  nonnegative  self-adjoint  solution  of  the  algebraic  Riccati  equation 

A*T[  +  UA-  UBR-^Bm  +  Q  =  0.  (3.1) 

Here  r  is  the  unique  r-periodic  solution  of 

r{t)  +  {A*  -  UBR-^Byit)  -  =  0  (3.2) 

and  the  optimal  trajectory  Z  is  the  solution  of 

Z{t)  =  {A-  BR-^B*U)Z{t)  +  BR-^B*r{t)  -}-  J^{t)  . 

These  equations  (in  particular  (3.1),  (3.2))  are  infinite  dimensional  (i.e.,  in  7f)  and  hence 
aj'proximation  techniques  a:e  required  to  obtain  approximate  feedback  gains.  Using  a  stan¬ 
dard  Galerkin  approach,  one  typically  chooses  a  sequence  of  finite  dimensional  subspaces 
C  H  with  projections  and  defines  an  approximating  problem  in  of 

minimizing 

J''(u)  =  ij('  +  {Ru(t)M‘))v}dt 

subject  to  an  approximating  system 

Z^{t)  =  A^Z^{t)  +  B^u{t)  +  J^^{t) 

Z^{0)  =  Z^{t}  =  V^Z{0). 

The  solutions  are  given  by 

u^{t)  =  -R-^B^’[U^Z^{t)  -  r^(t)] 

Z^it)  =  (>l^  -  B^R-^B’^*n'^)Z'^{t)  -I-  B^R-^B^>^[t)  -f  r^{t) 
where  H^  is  the  unique  nonnegative  self-adjoint  solution  of 


and  is  the  unique  r-periodic  solution  of 

r^{t)  -b  -  n'^B^R-^B^')r^{t)  -  =  0. 

In  order  to  guarantee  the  convergence  Yl^V^Z  -*  UZ  for  Z  €  H,  r^{t)  —*  r{t),  and 
hence  the  convergence  of  u^{t)  to  u{i),  it  is  sufficient  to  impose  various  conditions  on  the 
original  and  approximation  systems.  These  hypotheuos  include  convergence  requirements  for 
the  uncontrolled  problem  as  well  as  the  requirement  that  the  approximation  systems  preserve 
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stabilizability  and  detectability  margins  uniformly.  A  fully  developed  theory  (see  [4])  is 
available  for  the  case  that  =  0  (in  this  case  the  tracking  variable  r  does  not  appear  in  the 
solution)  even  in  the  case  that  B  is  unbounded  in  the  sense  formulated  in  Section  1.  The 
theory  in  [4]  requires  rather  strong  damping  assumptions  on  the  second  order  system  (2.5) 
in  order  to  be  applicable.  Under  appropriate  cissumptions,  the  techniques  and  ideas  of  [3] 
and  [4]  can  be  used  to  treat  the  case  for  .F  ^  0  in  both  identification  and  feedback  control 
problems. 


4  Finite  Dimensional  Approximation 

An  advantageous  feature  of  the  state  space  approach  for  feedback  control  is  that  the  optimal 
control  can  be  implemented  using  various  approximation  techniques.  To  illustrate  the  ideas 
involved,  let  denote  the  1-D  basis  functions  which  are  used  to  discretize  the  beam 

and  let  m  =  (m*  +  1)  •  {ruy  +  1)  —  1,  denote  the  2-D  basis  functions  which  are 

used  in  the  cavity.  The  n  —  1  and  m  dimensional  approximating  subspaces  are  then  taken 
to  be  =  span and  =  span ,  respectively.  Defining  N  =  m-{-n  —  l, 
the  approximating  state  space  is  x  if”  and  the  product  space  for  the  first  order 

system  is  =  II ^  x  .  The  finite-dimensional  approximation  is  then  determined  by 
restricting  cr  to  x  where  a  is  given  in  (2.6).  This  yields  the  operator  —*■ 

where 

-A’S 

and  and  are  obtained  by  restricting  ai  and  (72  to  x  II^ .  We  observe  that  the 
restriction  of  the  infinite  dimensional  system  (2.5)  to  the  space  x  yields  for  'P  =  (^,  rj) 

=  /  ■(()(//, i-W, -2)  flV7+  /  fvd-l- 

''To  t  ‘'To 


1=1 


When  ^  is  chosen  in  IJ^  and  the  approximate  beam  and  cavity  solutions  are  taken  to  be 


n-l 


i=I 


and 


1=1 
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respectively,  this  yields  the  system 


MV{t)  =  AV{t)  +  +  F^{t) 

M^y^(O)  =  Vq 


where 


(0  = 


Here  denotes  the  N  x  1  = 

jP  • 

(m  +  n  —  1)  X  1  approximate  state  vector  coefficients  while  u{t)  =  (ui(t),  •  •  •  «*(t))  contains 
the  s  control  variables.  The  full  system  has  the  form 


0  ]  [ 


0 

0 


i9^(0) 


0  ]  r 

-Af  -A^ 


0 1  r  0 

u(t)  + 

''  f«{i) 


M''=diaglM.«  A/a, 

A/i'=cliag(A/,';,A/a, 


0  <1 

■^32  ■^22 


= 


-  4': 


f  "(<)  = 
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The  component  matrices  are  given  by 

KL,*  =  /„  =  X.  D-‘Bre‘B;d-r , 

Klu  =  L  .  [M&l  .  =  nB:B;d-, , 

[<■]„ = ’  Kl,.i = X. ■ 

=  -l/iBrBTdi ,  . 

K],.  =  l^CBID^BtD^B;d^ , 


Moreover,  the  vectors  ^  and  elements 


[Sn],  =  X  =  X.  ' 

[9»l , = X  ’  [^“1. = X.  • 


In  all  cases,  the  index  ranges  are  k,£  =  1, •  •  •  ,m  and  i,p  =  1, •  •  •  ,n  —  1.  The  patch  index 
j  ranges  from  1  to  s.  It  should  be  noted  that  the  matrices  and  are  symmetric  and 
positive  definite  by  construction.  The  matrix  A2  has  a  symmetric  block  and  a  skewsymmetric 
block  and  the  eigenvalues  of  A2  are  real  and  nonnNigative. 

With  the  bzises  and  chosen,  the  finite  dimensional  theory  outlined  in 

the  Icist  section  holds  with  the  various  finite  dimensional  operators  replaced  by  appropriate 
matrices.  Specifically,  the  finite  dimensional  control  problem  is  then  to  find  u  6  Z/^(0,r) 
which  minimizes 

{{QV(‘),y''m)K''  +{Ru{t)Mt))R-}dt  ,  iV  =  m+n-l 
where  is  nonnegative  definite  and  t/^  solves 

y^{t)  =  /I  V(0  +  +  F^{i) 

s"(0)  = 

Here  A’^  =  =  (iV/'v)"*  5^^,  F^(0  =  with  the  initial 

condition  ijq  =  Po  .  The  optimal  control  is 

u^{t)  =  R-^  {B^f  [r'^CO  -  n^y^(0]  (4.2) 
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where  11^  is  the  solution  to  the  algebraic  Riccati  equation 

+  <5^  =  0.  (4.3) 

Since  denotes  the  matrix  representation  for  the  operator  Q^,  a  suitable  choice  for  is 

0 


=  V 


where  the  diagonal  matrix  V  is  given  by 

V  =  diag  [dir,d2r-\d3i^,d4P-^] 


Here  7** ,  =  m,  n  —  1 ,  denotes  a.kxk  identity  and  the  parameters  d,-  are  chosen  to  enhance 
stability  and  performance  of  the  feedback.  The  sxs  diagonal  matrix  R  contains  the  positive 
control  weights  and  has  entries  r,-,,z  =  For  the  regulator  problem  with  periodic 

forcing  function  F^{t),  r^{t)  solves  the  linear  differential  equation 


r^{t)  =  -  [a^  -  B^R-^{B^fU^Y  r^{t)  + 
r^(C)  =  r^(T) 

while  the  optimal  trajectory  is  the  solution  to  the  linear  differential  equation 

y^it)  =  [a^  -  y^{t)  +  B^R-^B^fr^{t)  +  F^{t) 

2/^(0)  =  2/^(r)  . 


(4.4) 


(4.5) 


5  Specific  Approximations  and  Numerical  Results 

We  next  turn  to  a  discussion  of  specific  choices  of  basis  functions  in  the  general  formulation  of 
the  approximation  schemes  in  the  Icist  section.  We  shall  also  present  numerical  results  from 
related  computations.  When  choosing  bases  for  the  finite  dimensional  subspaces  and 
in  a  control  setting,  one  must  weigh  criteria  such  cis  smoothness  requirements,  uniform 
preservation  of  exponential  stability  of  approximating  systems  (see  [5]),  accuracy,  sparsity 
of  system  matrices  and  ease  of  implementation. 

From  energy  considerations,  it  follows  that  the  system  (2.3)  is  dissipative;  hence  all  the 
eigenvalues  lie  in  the  left  half  plane.  The  model  of  Section  1  includes  no  medium  damping 
however,  and  hence  the  energy  dissipation  in  the  cavity  results  exclusively  from  the  boundary 
(Kelvin- Voigt  damping  in  the  beam)  thus  making  the  system  (2.3)  only  weakly  damped. 
In  spite  of  the  lack  of  strong  damping,  numerical  tests  have  indicated  that  when  physically 
relevant  parameters  are  used  in  the  model,  the  system  (2.3)  is  exponentially  stable  (the  fixed- 
end  boundary  conditions  on  the  beam  make  difficult  a  thorough  analytical  analysis  of  the 
eigenstructure).  When  considering  various  methods  of  discretizing  the  problem,  one  would 
like  to  choose  schemes  which  uniformly  preserve  the  exponential  decay  rate  as  the  dimension 
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of  the  approximate  system  (4.1)  increases.  This  can  be  easily  checked  by  determining  whether 
or  not  there  exists  a  uniform  margin  for  increasing  N  between  the  open  loop  eigenvalues  of 
the  system  matrix  ^4^  in  (4.1)  and  the  imaginary  axis. 

When  considering  the  control  problem,  one  is  also  concerned  with  the  preservation  of 
uniform  stabilizability  and  detectability  margin.^  for  the  closed  loop  approximation  systems. 
Hence  care  must  also  be  taken  so  that  approximation  schemes  are  chosen  so  as  to  preserve 
a  uniform  margin  between  the  closed  loop  eigenvalues  of  and  the 

imaginary  axis.  Numerical  schemes  which  satisfy  these  various  criteria  will  now  be  discussed. 

Cubic  splines  were  used  as  a  basis  for  since  they  satisfy  the  smoothness  requirement 
as  well  as  being  easily  implemented  when  adapting  to  the  fixed-end  boundary  conditions 
and  patch  discretizations.  For  a  given  positive  integer  n,  a  uniform  partition  Wcis  taken  with 
the  gridpoints  a:"  =  ^a,  i  =  0, l,-*-,n.  If  ^  is  used  to  denote  the  standard  cubic 

spline  basis  corresponding  to  this  partition  (see  [18],  page  79),  then  the  basis  functions  for 
the  beam  discretization  were  taken  to  be 

5"  =  B^-  2B^  -  23!!.^ 

Bf  =  Bf  ;  e  =  2, 3,  •  •  • ,  n  —  2 

It  is  readily  seen  that  these  basis  functions  satisfy  the  essential  boundary  conditions;  that  is, 


5^(0)  =  DB^iO)  =  (a)  =  DB^{a)  =  0 


for  i  =  1,2, •  •  •  ,n  —  1.  As  mentioned  previously,  the  corresponding  n  —  1  dimensional 
approximating  subspace  is  then  given  by  =  span  and  the  approximate  beam 

solution  is  taken  to  be 

f=i 

With  this  choice  of  basis  functions,  the  matrices  and  are  easily  constructed  and 
are  7-banded.  It  should  be  noted  that  a  Tau-Legendre  discretization  was  also  considered 
for  the  beam  but  had  the  disadvantage  of  the  loss  of  four  equations  due  to  the  constraints 
mandated  by  the  fixed-end  boundary  conditions  (see  [13]  for  a  discussion  of  Tau  methods). 

The  bases  that  were  considered  for  the  cavity  discretization  included  tensored  one¬ 
dimensional  Legendre  polynomials,  tensored  linear  splines  and  finite  elements.  The  methods 
of  system  formulation  as  well  as  the  advantages  and  disadvantages  of  each  can  be  summa¬ 
rized  as  follows.  Consider  first  the  Legendre  basis.  Let  P“(x)  and  P/(y)  denote  the  standard 
Legendre  polynomials  that  have  been  scaled  by  transformation  to  the  intervals  [0,a]  and 
[0,^],  respectively.  The  basis  functions  {P,^}  for  the  cavity  are  then  defined  as 

BTj{^,y)  =  Pi{x)Pj{y)  for  i  =  0, 1,- •  •  ,mx,  j  =  0, 1, •  •  •  ,mj„  z -j- j  ^  0, 

where  m  =  {nix  + 1)  •  (wy  + 1)  —  1.  The  condition  z  ^  0  eliminates  the  constant  function 
thus  guaranteeing  that  the  set  of  functions  is  suitable  as  a  bcisis  for  the  quotient  space.  For 
definiteness,  the  basis  functions  are  ordered  by  assuming  that  z  varies  for  each  fixed  j  which 
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is  analogous  to  a  left  to  right,  bottom  to  top  ordering.  Notice  that  because  natural  boundary 
conditions  occur  on  all  sides  of  the  cavity,  one  does  not  have  to  employ  a  Tau  method;  that 
is,  the  method  is  simply  a  Galerkin  scheme  without  modification  of  the  basis  elements  to 
satisfy  some  essential  boundary  conditions. 

The  component  matrices  and  can  then  be  succinctly  described  as  follows.  Let 
the  fundamental  (mj,  +  1)  x  {rux  +  1)  matrices  and  be  defined  as 

WTh  =  [ 

jO 

[KTh  =  J‘ DPnx)DP;(x)dx 

with  similar  definitions  for  M^,  K^.  Using  the  tensor  properties  of  the  2-D  bzisis,  we  can 
form  the  matrices  and  defined  by 

=  Mp  0  0 

The  ordering  in  the  above  definition  depends  on  the  ordering  of  the  basis  functions.  The 
matrices  and  are  obtained  by  removing  the  first  row  and  first  column  of  and 
to  reflect  the  deletion  of  the  constant  function  from  the  basis  set.  Note  that  with 
this  definition,  both  matrices  are  very  easily  constructed  and  that  the  mass  matrix  is 
diagonal;  hence  the  inverse  is  trivial  to  calculate.  Although  the  matrix  is  not  sparse,  it 
has  a  well-defined  structure  due  to  its  tensor  product  nature  and  the  fact  that  and  M/" 
are  diagonal.  It  too  can  be  efficiently  inverted  when  one  takes  advantage  of  this  structure.  In 
the  ceise  that  pj  is  constant,  the  stiffness  matrix  Afi  can  be  constructed  in  the  same  manner 
as  Mil  ^he  tensor  product  structure  can  be  used  advantageously  both  when  solving  the 
Riccati  equation  (4.3)  and  the  ODE  systems  (4.4)  and  (4.5). 

In  order  to  use  a  tensored  linear  spline  or  finite  element  basis  in  the  cavity,  some  constraint 
must  be  applied  in  order  to  guarantee  that  is  a  quotient  space  (one  cannot  simply  drop 
the  constant  function  as  was  done  with  the  tensored  Legendre  polynomials).  One  such 
constraint  which  is  commonly  used  is  the  requirement  that 

/  (j)^ {t,x^y)(kij  =  0  . 

If  denote  the  standard  tensor  product  linear  spline  bcisis  (see  page  129 

of  [18]),  then  the  integral  constraint  leads  to  the  quotient  basis  where 

=  Br{x,y)  -  4afBi”‘(x,y) 

with  a<  =  depending  upon  whether  the  function  is  a  corner  basis  function,  a 

side  basis  function  or  an  interior  basis  function,  respectively.  As  a  result  of  the  modifica¬ 
tions  needed  to  obtain  a  quotient  space  bcisis  with  the  tensored  linear  splines,  the  matrices 
Mil,  M^2  ^11  hence  one  loses  the  structural  advantages  obtained  with  the 

Legendre  basis. 
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Similar  modifications  must  be  made  when  using  a  finite  element  basis  in  a  quotient  space 
with  the  result  that  the  system  matrices  are  also  full  in  that  case.  Moreover,  the  lower 
order  accuracy  of  the  splines  and  finite  elements  necessitates  the  use  of  a  larger  number  of 
basis  functions  and  hence  larger  matrices  in  order  to  match  the  accuracy  of  the  Legendre 
polynomials.  The  fact  that  the  Legendre  basis  yields  smaller,  structured  matrices  than 
those  obtained  with  the  linear  splines  and  finite  elements  is  important  but  not  crucial  in  the 
problem  under  consideration  since  the  cavity  is  only  two  dimensional  and  hence  matrix  sizes 
are  reasonably  small.  This  issue  will  become  much  more  critical  when  considering  the  3-D 
problem  of  interest  because  of  the  large  matrix  sizes  which  will  be  encountered. 

As  discussed  earlier,  a  final  item  which  should  be  considered  when  choosing  a  means 
of  discretizing  the  control  problem  is  whether  or  not  the  approximation  scheme  effects  a 
uniform  preservation  of  exponential  stability  for  the  open  and  closed  loop  approximating 
systems.  This  issue  is  illustrated  by  the  results  in  the  Example  5.1. 

The  problem  under  consideration  in  Examples  5.1  and  5.2  is 


4)tt  =  c^A<^  (x, j/)  €  D  >  0  , 

V^-n  =  0  (x,i/)  €  r  >  0  , 

^(t,x,0)  =  -Wt{t^x)  0<x<.6,f>0, 
oy 


dx^ 
dx"^ 

\ 


EI—u{t)  [H{x  -  a„)  -  H{x  -  a, 2)] 


0  <  X  <  .6  ,  i  >  0  , 


0)  =  0)  =  w{t,  .6)  =  -^(t,  .6)  =  0  t>  0  , 


(5.1) 


^(0,x,y)  =  </>t(0,x,y)  =  t/;(0,x)  =  wt(0,x)  =  0 


where 

f{t,x)  =  2.04sin(1507r<)  . 

The  parameter  choices  a  =  .6  m,  £  =  1  m,  Pf  =  1.21  kgfmP,  =  117649  rn^lsec?^ 
pb  =  1.35  kg/m,  El  =  73.96  Nm'^,  cq!  =  .001  kg  m^/ sec,  =  82.9629,  T  =  .0005  m, 
^  =  1.9  X  10"^°  «ii  =  -25  and  or, 2  =  .35  are  physically  reasonable  for  a  .6  m  by  1  m 

cavity  in  which  the  bounding  end  beam  has  a  centered  piezoceramic  patch  covering  1/6  of 
its  length  (see  Figure  3).  The  beam  is  assumed  to  have  width  and  thickness  .1  m  and  .005  m, 
respectively.  The  quadratic  cost  functional  parameters  were  taken  to  be  f/i  =  (£2  =  (U  =  1, 
^3  =  10^  and  R  =  10~®  with  c/3  of  much  larger  magnitude  than  di,d2  or  d.^  to  emphasize 
the  penalization  of  large  pressure  variations.  Note  that  because  there  is  only  one  patch,  the 
control  weight  R  is  simply  a  positive  scalar. 

For  a  beam  with  the  above  dimensions  and  density,  the  natural  frequency  of  the  first 
mode  is  73.21  hertz  and  the  frequency  of  the  forcing  function  was  chosen  so  as  to  be  close 
to  this  value.  To  obtain  the  magnitude  2.04,  it  was  assumed  that  the  forcing  function  was 
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the  result  of  an  exterior  plane  wave  with  a  sound  pressure  level  of  120  dB  (which  forces  an 
interior  sound  pressure  level  of  98  dB). 


0  .6 


Figure  3.  Example  acoustic  chamber  with  one  piezoceramic  patch. 


Example  5.1 

In  this  example,  the  uniform  preservation  of  exponential  stability  for  the  open  and  closed 
loop  approximating  systems  is  examined.  For  n  =  =  m,,  =  5, 6, 7  and  8,  the  margins  of 

stability  for  the  open  and  closed  loop  systems  obtained  with  tensored  Legendre  polynomials 
and  tensored  linear  splines  are  listed  in  Tables  1  and  2,  respectively.  The  gains  needed 
for  the  closed  loop  system  were  calculated  via  Potter’s  method  (see  [17]).  For  each  n,  the 
locations  of  the  open  and  closed  loop  eigenvalues  obtained  with  the  Legendre  polynomials 
are  displayed  in  figures  4  and  5,  respectively.  When  plotting  the  eigenvalues  of  and 

—  B^R~^  \^j  those  eigenvalues  having  real  parts  with  magnitude  greater  than 
1  have  been  excluded  in  order  to  better  see  the  distribution  near  the  imaginary  axis.  Note 
that  a  uniform  margin  of  stability  is  maintained  between  both  the  open  and  closed  loop 
eigenvalues  and  the  imaginary  axis  for  both  sets  of  bases.  Results  similar  to  those  obtained 
with  the  Legendre  basis  were  obtained  when  finite  elements  were  used  as  a  basis  for 


Table  1.  Margin  betv/een  the  open  and  closed  loop  eigenvalues  and  the  imaginary  axis  with 
tensored  Legendre  polynomials. 


rux  =  rriy  n  max  {i?eA,  A  6  max  ^Re\,  X  e  a  -  B^R-^  } 


7  7  -.0200  -.0200 

8  8  -.0158  -.0290 


Table  2.  Margin  between  the  open  and  closed  loop  eigenvalues  and  the  imaginary  axis  with 
tensored  linear  splines. 


Figure  4.  Eigenvalues  of  for  n  =  rrix  =  ruy  =  5,6,7  and  8  with  tensored  Legendre 
polynomials. 


7  ' 

Figure  5.  Eigenvalues  of  —  B^R~^  n  =  7nj:  =  my  =  5, 6, 7  and  8  with 

tensored  Legendre  polynomials. 


As  seen  in  Example  5.1,  a  larger  margin  of  stability  is  maintained  in  both  the  open  and 
closed  loop  systems  with  the  linear  spline  basis  than  with  the  Legendre  basis;  hence  one  might 
conclude  that  the  linear  splines  are  the  basis  of  choice  when  solving  the  control  problem.  As 
noted  earlier  however,  one  must  also  weigh  factors  such  as  system  size,  accuracy  and  efficiency 
when  choosing  a  numerical  method.  Numerical  tests  have  indicated  that  in  spite  of  the  larger 
eigenvalue  margins  of  the  linear  splines,  their  performance  when  used  in  the  control  problem 
is  nearly  identical  to  that  obtained  with  the  Legendre  polynomials.  Moreover,  because  of  the 
lower  order  accuracy  of  the  splines,  a  larger  number  of  basis  functions  is  needed  to  obtain 
suitable  accuracy  thus  leading  to  matrix  dimensions  that  are  almost  twice  those  resulting 
from  the  Legendre  discretization.  Finally,  as  noted  earlier  in  this  section,  the  matrices 
obtained  with  the  Legendre  discretization  are  much  more  structured  than  those  obtained 
with  finite  elements  or  linear  splines  hence  making  Legendre  implementation  more  efficient 
than  the  other  cases.  Results  for  the  LQR  control  problem  for  (5.1)  with  the  tensored 
Legendre  basis  for  the  cavity  and  cubic  splines  for  the  beam  with  iiix  =  niy  =  4  and  n  =  8 
are  reported  in  Example  5.2. 
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Example  5.2 

In  this  example,  the  effect  of  the  feedback  ct  ntrol  on  the  problem  for  system  (5.1)  is  de¬ 
scribed.  In  order  to  solve  for  the  optimal  control  and  trajectory,  it  is  necessary  to  solve  both 
the  trajectory  equation  (4.5)  and  the  tracking  equation  (4.4).  Because  numerical  evidence 
indicated  that  both  unconstrained  solutions  were  roughly  periodic  with  period  r  =  1/75,  the 
problems  were  solved  ais  initial  value  problems  with  starting  values  j/(0)  =  0  and  r(10/75)  =  0 
rather  than  as  free  boundary  value  problems.  The  choice  for  initial  state  is  physically  rea¬ 
sonable  while  the  choice  to  integrate  backwards  in  time  in  (4.4)  is  made  to  reduce  numerical 
instability  when  solving  the  ODE  system  for 

The  uncontrolled  and  controlled  approximate  acoustic  pressures  {p^  =  Pj<f>^)  at  the 
point  (X,Y)  =  (.3,  .1)  are  plotted  in  Figure  6  for  the  time  interval  [0, 10/75].  Similar  plots 
for  the  approximate  beam  displacement  at  X  =  .3  are  given  in  Figure  7.  The  uncontrored 
solutions  exhibit  a  beat  phenomenon  which  results  from  the  fact  that  the  frequency  of  the 
forcing  function  is  slightly  .ater  than  the  natural  frequency  of  the  first  mode  of  the  beam. 
After  a  transient  interval,  the  controlled  rolutions  are  periodic  and  are  maintained  at  a 
level  which  is  approximately  10%  of  that  found  in  the  uncont!  olled  case  (note  the  scales  in 
Figures  6  and  7).  This  produces  an  interio.  sound  pressure  level  of  77  dB  which  is  a  21  dB 
reduction.  To  further  illustrate  the  state  reduction  with  feedback  control,  the  uncontrolled 
and  controlled  acoustic  pressures  at  the  times  T  =  1/75,  2/75,  6/75  and  10/75  are  plotted 
in  Figures  8-11,  respectively.  The  two  dimensional  plots  in  each  figure  show  spatial  slices 
of  the  uncontrolled  and  controlled  pressures  at  A"  =•  .3,  0  <  ?/  <  1.  Figures  12  and  13 
contain  plots  of  the  uncontrolled  and  controlled  beam  displacements  at  the  times  T  =  6/75 
and  T  =  10/75,  respectively.  The  results  in  Figures  8  -  13  are  reoresentative  of  those  found 
throughout  the  time  interval  (0, 10/75]  and  in  conjunction  with  Figures  6  and  7,  demonstrate 
that  the  pressure  and  beam  displacement  are  uniformly  reduced  and  maintained  at  a  very 
low  level  of  magnitude  in  spite  of  the  periodic  forcing  function. 

The  controlling  voltage  u{t)  is  plotted  in  Figure  14.  As  expected,  it  is  periodic  with 
period  1/75.  It  should  be  noted  that  the  magnitude  of  u(t)  remains  less  than  60F  which  is 
a  physically  reasonable  voltage  to  put  into  the  piezoceramic  patches. 

As  mentioned  in  the  last  section,  the  choice  of  the  quadratic  cost  functional  parameters 
di  —  di  and  R  influences  the  control  stability  and  performance  of  the  feedback.  In  this 
problem,  the  emphasis  is  on  the  the  reduction  of  variations  in  the  acoustic  pressure;  hence 
dz  was  taken  to  be  larger  than  c/i,  dz  or  d^.  It  should  be  noted  that  this  choice  of  parameters 
does  not  exclude  the  control  of  the  other  state  variables;  in  fact,  the  beam  displacement  is 
significantly  reduced  as  seen  in  Figures  7,  12  and  13.  Since  the  parameter  i?  is  a  penalty 
term  for  u(t),  more  control  of  the  state  variables  can  be  effected  by  choosing  R  smaller.  The 
tradeoff,  however,  is  an  increase  in  the  voltage.  Hence  one  must  weigh  the  amount  of  state 
reduction  desired  against  the  amount  of  voltage  which  can  be  put  into  the  patches. 

The  amount  of  control  is  also  directly  influenced  by  patch  size,  placement  and  the  number 
of  patches  being  used.  In  the  examples  that  we  have  observed,  the  best  results  were  obtained 
with  one  centered  patch,  and  we  have  noticed  that  the  amount  of  control  obtained  increases 
with  increasing  patch  length.  Thus  one  must  weigh  the  amount  of  control  desired  against 
physical  limitations  on  the  size  of  the  piezoceramic  patches  being  used. 
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Figure  6.  Uncontrolled  and  controlled  pressures  at  the  point  (X,y)  =  (-3,  .1)  throughout 
the  time  interval  [0, 10/75]. 


Figure  7.  Uncontrolled  and  controlled  beam  displacements  at  the  point  X  -  .3  throughout 
the  time  interval  [0, 10/75). 
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Figure  12.  Uncontrolled  and  controlled  beam  displacements  at  T  =  6/75 
- uncontrolled  displacement, - controlled  displacement. 

xlO-*  T=  10/75 
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Figure  13.  Uncontrolled  and  controlled  beam  displacements  at  T  =  10/75 
- uncontrolled  displacement, - controlled  displacement. 
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6  Conclusion 


For  the  2-D  acoustic  problem  involving  the  transmission  of  exterior  noise  into  an  interior 
cavity  via  fluid/structure  interactions,  a  model  set  of  differential  equations  has  been  devel¬ 
oped.  Control  is  implemented  in  the  model  via  piezoceramic  patches  on  the  beam  which  are 
excited  in  a  manner  so  as  to  produce  pure  bending  moments.  By  writing  the  resulting  system 
as  an  abstract  Cauchy  equation,  the  problem  of  reducing  interior  pressure  fluctuations  can 
be  posed  in  the  context  of  an  LQR  time  domain  state  space  formulation  and  approximation 
schemes  which  are  suitable  for  this  theory  are  presented. 

For  one  typical  patch  configuration,  examples  are  given  which  demonstrate  the  stabiliz- 
ability  of  the  open  and  closed  loop  systems  under  approximation  as  well  as  the  reduction  of 
cavity  pressure  and  beam  displacement  when  the  feedback  control  is  invoked.  The  examples 
show  that  input  of  the  optimally  controlling  voltage  u{t)  uniformly  reduces  both  the  pressure 
and  the  beam  displacement  and  maintains  them  at  a  very  low  level  of  magnitude  throughout 
the  time  interval  of  interest. 

As  mentioned  in  the  example  section,  the  amount  of  control  obtained  is  directly  influenced 
by  patch  size,  placement  and  the  number  of  patches  being  used.  Initial  results  have  indicated 
that  for  a  uniform  periodic  forcing  function,  the  best  results  can  be  obtained  with  one 
centered  patch  with  the  amount  of  control  increasing  with  increasing  patch  length.  Further 
computational  studies  are  currently  being  conducted  to  determine  the  effect  of  various  pat^'h 
configurations  and  forcing  functions  on  the  decibel  reduction  in  the  cavity. 
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